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Abstract 

We consider solutions to the linear wave equation O g <f> — on a non- 
extremal maximally extended Schwarzschild-de Sitter spacetime arising 
from arbitrary smooth initial data prescribed on an arbitrary Cauchy 
hypersurface. (In particular, no symmetry is assumed on initial data, 
and the support of the solutions may contain the sphere of bifurcation of 
the black/white hole horizons and the cosmological horizons.) We prove 
that in the region bounded by a set of black/ white hole horizons and 
cosmological horizons, solutions (j> converge pointwise to a constant faster 
than any given polynomial rate, where the decay is measured with respect 
to natural future-directed advanced and retarded time coordinates. We 
also give such uniform decay bounds for the energy associated to the 
Killing field as well as for the energy measured by local observers crossing 
the event horizon. The results in particular include decay rates along the 
horizons themselves. Finally, we discuss the relation of these results to 
previous heuristic analysis of Price and Brady et al. 
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1 Introduction 

The introduction of a positive cosmological constant in the Einstein equations 
of general relativity gives rise to a wide variety of new interesting solution space- 
times, in particular, spacetimes containing both "black hole" and "cosmological" 
regions. As in the case of black- hole spacetimes with vanishing cosmological con- 
stant, the stability of these spacetimes as solutions to the Einstein equations is a 
fundamental open problem of gravitational physics. Yet even the simplest ques- 
tions concerning the behaviour of linear waves on such spacetime backgrounds 
today remain unanswered. In this paper, we initiate in the above context the 
mathematical study of decay for solutions to the linear wave equation. 

The simplest family of black-hole spacetimes with positive cosmological con- 
stant is the so-called Schwarzschild-de Sitter family. If the cosmological constant 
A > is considered fixed, this is a 1-parameter family of solutions (A4,g) to the 
Einstein vacuum equations 

R»v ~ T.9livR = -h-9p.v, (1) 
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with parameter M, called the mass. We shall consider only the non-extremal 
black-hole case, corresponding to parameter values 



< M < 



1 



3VA' 



(2) 



As with the Schwarzschild family, the first manifestation of the Schwarzschild- 
de Sitter family of solutions was an expression for the metric in local coordi- 
nates, in this case first published in 1918 by Kottler [17], and independently by 
Weyl [20], in the form 



- 1- 



2M 



1 - 



2M 



- rAr 2 



dr 2 



r 2 da R 



(3) 



Here dcr§2 denotes the standard metric on the unit 2-sphere. The global struc- 
ture of maximal spherically symmetric vacuum extensions of such metrics was 
only understood much later [5] [TS] [13] based on the methods of formal Pen- 
rose diagrams introduced by B. Carter. In fact, maximally extended spherically 
symmetric vacuum spacetimes {M , g) with various different topologies can be 
constructed, all of which equally well merit the name "Schwarzschild-de Sitter 
with parameter M and cosmological constant A". Such solutions (M,g) all 
share the property that the universal cover Q of the 2-dimensional Lorentzian 
quotient Q = M/SO(3) consists of an infinite chain of regions as depicted in 
the Penrose diagram below: 



The results of this paper do not depend on the topology, and for definiteness, 
one may assume in what follows that the name "Schwarzschild de-Sitter" and 
the notation {M. , g) refer to the spacetime with quotient precisely the universal 
cover depicted above. 

It is then the wave equation 

u g 4> = (4) 

on this background (M.,g) whose mathematical study we wish to initiate here. 
There is already a rich body of heuristic work on this problem in the physics 
literature. (See Section [L4"l bclow for a discussion.) The motivation for the study 
of ([4} in the present context is multifold. In particular, as in the case of vanishing 
cosmological constant, studied in our previous [12] , we believe that proving 
bounds on decay rates for solutions to ([4]) is a first step to a mathematical 
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understanding of non-linear stability problems for spacetimes containing black 
holes, that is to say, to the problem of stability in the context of the dynamics 
of ([1]). For a more detailed discussion, we refer the reader to the introductory 
remarks of [12] , 

1.1 The initial value problem for the wave equation 

We are interested in solutions of ((4|) arising from suitably regular initial data 
prescribed on a Cauchy surface E of M. . For future applications to non-linear 
stability problems, it is crucial that all assumptions have a natural geometric 
interpretation independent of special coordinate systems. Moreover, our pri- 
mary concern in this paper is the region T> bounded by a set of black/white hole 
horizons 7i + U Or and cosmological horizons Tt U TC : 

V = clos(j-(H + un + )n j + (ruH")) (5) 

as depicted 

belo^ 



P\~~ r = 












a--"" r = ~ 




r = oo 



By causality, the global behaviour of <p in D can be understood independently 
of the behaviour near r = and r = oo. The behaviour in say T> n J + (E) 
is completely determined by the behaviour of appropriate initial data on E D 
review briefly in the next paragraph the solvability and domain of 
dependence property for the initial value problem for Q. 

Let E C M. be a smooth Cauchy surface and let n 11 denote the future- 
directed unit normal of E. For s > 1, let <p be an Hf oc (E) function and ip : E — > K 
an i/ 1 s o ~ 1 (E) function. Then there exists a unique global solution 4> : M. — > M 
of n g (j) = such that for all smooth spacelike hypersurfaccs E with future 
directed unit normal h, G H\oc ("- < / , )|s / € -^loc ' and <f>\z = tp, ncf>\s = <p- 
Moreover, if K C A4 is closed and cj>i, (f>2 are two such solutions corresponding 
to data (<pi,(pi), (<P2,^2) such that pi\-E>nK = ^Is'nK, V>i\s'nK = ^Is'rw, 
then 0i = 02 on7W\(J+(E\ii:)UJ"(E\iY)). In particular, setting K = J~{T>), 
we obtain that <j>i = 4>2 on J + (E) n 2?. 

We employ in this paper the standard notation of Lorentzian geometry (e.g. J + , J - , 
etc.), and Penrose diagrams. See |15| . 

2 Note, as depicted, that S n J~ (T>) is not necessarily £ n T>. 
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1.2 The main theorem 

1.2.1 Norms on initial data 

Since our results will be quantitative, we need to introduce relevant norms on 
the compact manifold with boundary En J - (T>). Let || • || denote the Ricmannian 
L 2 norm on E n J~(T>). This induces a norm on sections of the tangent bundle, 
a norm we will denote also by || • ||. If if £ Hl oc (E), if £ L 2 oc (S), then let us 
denote by <f> the unique solution of d g <fi = corresponding to initial data {if, if). 
Let us define now for all real s > the quantity 

E.(<p,<j>) = ||V S ^|| 2 + ||0|| 2 +J2 r2se2s \\V^i\\ 2 +r 2s e 2s \\nfa\\ 2 , (6) 

where 4>g denotes the projection of cf> to the fth eigenspace of //\, i.e. the £'th 
spherical harmonic of </>. The function r is discussed in Section [2j If itself 
is spherically symmetric, then we may replace 4>i, ncf>i be ifi and ifj>, and the 
above expression is a sum of integrals on initial data. For general E, a sufficient 
condition for the finiteness of © is that if £ iJ^+^E), if £ iff oc (E). 

In the case m > an integer, we can characterize E m geometrically as 
follows. Let fij, i = 1,...,3 denote a basis of Killing fields generating the 
Lie algebra so(3) associated to the spherical symmetry of {A4 , g) . We call fli 
angular momentum operators. It easily follows that 

v m {f,f)~ y, E ii v ^ 

jpi)---Pm— i=0,l l<ii,...i m _i<3 

+ KO£ 

Again, if E itself is spherically symmetric, we may replace <p with f in the first 
term, and remove the n from the second, replacing <j> with ip. 

1.2.2 First statement of the theorem 

The main result of this paper is contained in the following 

Theorem 1.1. Let {M,g) denote the Schwarzs child- de Sitter spacetime with 
parameter M and cosmological constant A satisfying (0) and let E be a Cauchy 
surface for M.. Let T> C M. denote a region as defined in and let s > 0. 
Then, there exist constants C s depending only on s, M , A, and the geometry 
of E n J~(T>) such that for all solutions <f> of the wave equation n g = on 
M such that E s {ip,ip) is finite, where if = <f = ^0|s, and for all achronal 
hypersurfaces E' C T> fl J + (E'), the bound 

f T^Wr? <C s E 8 {ip,if){v + {T l ')- s + u + {T,')- s ) (7) 

holds, where u and v denote fixed Eddinqton-Finkelstein advanced and retarded 
coordvnateB, u + ^ max{u, 1}, u + {E) = inf u + {x), etc, denotes the 

3 See Section [2] Although these coordinates are only defined in 23°, the statements (0, JH} 
can be interpreted in all of T> in view of conventions JIQP — . 



•■<r»i 
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Killing field coinciding in the interior ofT> with J^, T^ u ((f>) denotes the standard 
energy-momentum tensor, and n v is the future- directed unit normal wherever £' 
is spacelike, in which case the integral is taken with measure the induced volume 
form^ 

In addition, |?J holds if is replaced by the vector field 7V M defined in 
Section® If s > 1, then the pointwise bound 

|^-^|<C s Ej(^ i j)(«^ 1 +^) (8) 

holds in J + (S) r\T>, where <p is a constant satisfying 

|0| < sup |0(x)| + Co ( V> #. 
— ies 

In particular, the theorem applies to arbitrary smooth initial data 

if G C°°(S), ip G C°°(S), where s can be taken arbitrarily large. There are 
no unphysical assumptions regarding vanishing of <p at the sphere of bifurcation 
of the horizons, i.e. at the sets H + P\H~ and H + C\H . The decay rates ©, (JSJ) 
are uniform, i.e. they hold up to and including the horizons, setting u + = oo or 
v + = oo. In particular, E' in ([9]) can be taken (as depicted below) 




r = oo 



to contain subsets of Tt + and/or TL . 

1.2.3 Comparison with the Schwarzschild case 

The statement of Theorem 11.11 should be compared with the results of our 
previous |10l 112] concerning the wave equation on a Schwarzschild exterior. 

Recall that in the region r > 2M, the Schwarzschild metric is given by the 
expression ([3]) for A = 0, M > 0. The Penrose diagramme of the closure of this 



4 A correct interpretation of and the measure of integration for general achronal S' can 
be derived by a limiting procedure. 
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region in the maximally extended Schwarzschild spacetime is given below: 




In |10j , an analogue of |7|) is proven for all s < 6 and spherically symmetric 
initial data. Modulo an e, this result is expected to be sharp, as it not expected 
to be true for s > 6, in view of heuristic arguments due to Price [19] . 

In |12j , an analogue of ([7]) is proven for s = 2 for arbitrary, not necessarily 
spherically symmetric, initial data. 

In view of the fact that solutions of the wave equation vanish on X + , the 
results of [T^] allow one to obtain the uniform pointwise decay rate \<f>\ < Cv 7 . 
As a uniform decay bound in v, this decay rate is in fact sharp. 



1.2.4 Second statement of the theorem 

The loss of angular derivatives in the result of Theorem ll . ll can be more precisely 
quantified by decomposing (f> into spherical harmonics. Each spherical harmonic 
4>i decays at least exponentially, but the bound on the exponential rate obtained 
here decreases inverse quadratically in the spherical harmonic number. We have 
the following 

Theorem 1.2. Let (A4,g), S, T> be as in Theorem \l.ll Then there exists a 
constant c depending only on M and A, and Co depending only on M , A, and 
the geometry o/Efl J~(T>), such that for all <^n solutions of the wave equation 
on A4 with spherical harmonic number £ with 

E (^, fa) = || V^|| 2 + ||<yi£|| 2 < oo 

and all achronal hypersurfaces £' C J + (S) n T>, the bound 

f T^Wn" < C Vo(w^i)(e- 2cv +^ )/e2 + e -WS'>/*) (9 ) 

holds for all £ > 0, and, again as before, also with T M above replaced with 
defined in Section^ In addition, the pointwise bounds 

|&(u,v)| < C E ^,^)(e-™^ 2 +e- cu +^ 2 ) 

for £ > 1 , and 

\Mu, v) -0| < Co Ei(<po, ^o)(e-™ +/£2 + e-™+/ £2 ), 
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for £ = 0, hold in J + (E) n T>, where <p is a constant satisfying 

\<t>\ < mf |0o(z)l+CoE o 4 (^ o ,¥>o). 
The above theorem can easily be seen to imply Theorem ll.il 

1.3 Overview of the proof 

In this paper, we insist on a framework of proof that in principle may have 
relevance to the non-linear stability problem, that is to say, the problem of the 
dynamics of U) starting from initial data close to those induced on a Cauchy 
hypersurface S of Schwarzschild-dc Sitter. This leads us to try to exploit com- 
patible currents. In this section, we will describe this general approach, and 
the natural relation of the currents we will define with various geometric and 
analytical aspects of the problem at hand. 

1.3.1 Vector fields and compatible currents 

For quasilinear hyperbolic systems (like ([I}) in 3 + 1 dimensions, all known 
techniques for studying the global dynamics are based on L 2 estimates. In the 
Lagrangian case, the origin of such estimates can be understood geometrically 
in terms of compatible currents (see Christodoulou [3])- These are 1-forms J M 
such that at each point both J M and the divergence K = V M J M depend 

only on the 1-jet of 4>. An important class of these are the currents obtained 
by contracting the energy momentum tensor T Ml/ with an arbitrary vector field 
V*. See Section [3] for a discussion in the context of the linear wave equations 
d g (j) = studied here. All estimates in this paper are obtained by exploiting 
the integral identities 

/ K = [ J„n" (10) 

corresponding to compatible currents of the form and straightforward mod- 
ifications thereof J M = + ■ ■ ■ , where the vector fields V are directly related 
to the geometry of the problem, and the region 1Z is suitably chosen. 

1.3.2 The photon sphere and the currents Jrf 

The timclike hypersurface r = 3M is known as the photon sphere. This has the 
ominous property of being spanned by null geodesies. If additional regularity is 
not imposed, then it is clear by a geometric optics approximation that solutions 
of the wave equation can concentrate their energy along such geodesies for ar- 
bitrary long times, and one can thus not achieve a quantitative bound for the 
rate of decay in terms of initial energy alone. In particular, |(7j) cannot hold for 
s > if E s is replaced by E . 

It is truly remarkable that this obstruction arising from geometrical optics 
is captured, and quantified, by a current associated to a vector field V 
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of the form f{r*)d r - for a well-chosen function /Jf| The story is not entirely 
straightforward, however. The desired current is in fact not precisely of the 
form jY , but a modification thereof, to be denoted J x,s , which is associated in 
a well defined way to a collection of vector fields Xg = fi(r*)d r *. The current 
is defined by summing over currents J Xi > 3 which act on individual spherical 
harmonics (f>g. 

The current J X ' 3 yields a nonnegative K X ' 3 , modulo an error term supported 
near the horizons. In a first approximation, we may pretend that in fact K X ' 3 > 
0, but degenerates (in regular coordinates) near the horizon. The identity (fTU|) 
can then be used as an estimate for its left hand side, in view of the fact that 
its right hand side will in fact be bounded by the flux of J T , for the Killing field 
T, which is conserved. The role of the photon sphere will be exemplified by the 
degeneration at r = 3M of the quantity controlled by this spacetime integral. 

In order to obtain decay results from the above, one would have to gain infor- 
mation about the quantity estimating the boundary terms-namely Jg-fiJ^n^, 
from the control of spacetime integral. The difficulty for this is that the space- 
time integral estimates one obtains degenerate at the photon sphere r = 3AI 
and at the horizons. This docs not allow one to control there. 

The problem at the photon sphere is cured by applying the estimate also to 
angular derivatives. It is here that the argument "loses" an angular derivative. 
It is this loss that leads to the form of decay proven in ([7|). 

The problem on the horizon, on the other hand, turns out to be illusionary. 
The horizon is in fact a very favourable place for estimating the solution! For 
this, we will need to consider the "local observer" vector fields Y, Y, to be 
described in the next section. 

1.3.3 The red-shift effect and the currents J„ and J Y 

The heuristic mechanism ensuring decay near the horizons has been understood 
for many years, and is known as the red shift effect. This is typically described 
in the language of geometric optics. If two observers A and B cross the event 
horizon at advanced times Va < Vb, and A sends a signal to B at a certain 
frequency, as he (A) measures it, then the frequency at which B receives it is 
exponentially damped in the quantity vb — va ■ 

It turns out that this exponential damping property can be captured by 
the integral identities ([TO)) corresponding to the currents J Y and J Y associated 
to vector fields Y, Y, defined in Section [8l These vector fields are supported 
near the horizons TC + , H , respectively. The estimates (|10|) corresponding to 
the currents J Y , J Y fulfill the double role of (a) correcting for the error region 
where K x ,3 < by dominating this term near the horizon by K Y + K Y and (b) 
controlling the spacetime integrated energy measured by local observers near 
the horizon. The choice of Y, Y is delicate, because there is an "error region" 
where K Y + K Y < 0, which must be controlled with the help of the currents of 

5 This insight, in the case of the wave equation on the Schwarzschild solution, is originally 
due to Blue and Soffer [2]. See, however, [3]. 
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the previous section. The use of the currents corresponding to X, Y and Y are 
thus strongly coupled. 

1.3.4 Comparison with the Schwarzschild case 

To see the above arguments in context, the reader may wish to compare with 
our previous [12) . where versions of the currents J x , J Y are also employed. 
The relation of our arguments with the physical mechanisms at play are in fact 
much clearer in the present paper, than in [12j . This is due on the one hand to 
the absence here of the Morawetz-type vector field (denoted K in [Hj), and, on 
the other hand, to the relative simplicity here in the construction of the current 
J a ■ We give here some comments on these points. 

The Morawetz vector field employed in [12] is a highly unnatural quantity 
at the horizon from the geometric point of view. On the other hand, in view 
of its weights, it somewhat magically captures a polynomial (as opposed to the 
proper exponential) version of the red shift. The pointwise decay rates achieved 
via K at the horizon are worse than the decay rates away from the horizon, but 
sufficient if one is only interested in the behaviour of the solution away from the 
horizon. (See also [4j.) In our [12] . uniform decay rates up to the horizon were 
indeed obtained with the help of J Y . But these estimates could be obtained 
a posteriori. From the point of view of the non-linear stability problem, this 
decoupling appears to be an exceptional feature. It is in this sense that the 
scheme proposed in the present paper is perhaps more naturally connected to 
the geometry of general black holes. 

The second point to be made here concerns the construction of Jjf . In [12] , 
positivity of the analogue of what we denote here K X 3 relied on an unmotivated 
recentring and rescaling of the derivatives of the functions ft which obscured 
perhaps the fundamental connection with the photon sphere. Here, this con- 
nection appears much more clear. Of course, this is at the expense of having 
to bound — K X ' 3 from K Y and K Y . This should in no way be thought of as a 
disadvantage. The red-shift effect has a lot to offer. It should be used and not 
obscured. 

1.4 Discussion 

As noted above, the study of the asymptotic behaviour of solutions to U g <f) = 
on both Schwarzschild and Schwarzschild-de Sitter backgrounds has a long 
tradition in the physics literature. In the Schwarzschild case, the pioneering 
heuristic study is due to Price [19]. See also [14]. For the Schwarzschild-de 
Sitter case, there is numerical work of Brady et al [B], the subsequent [7], and 
references therein. 

The above studies are based entirely upon decomposition of 4> into spherical 
harmonics. The results of these heuristics or numerics are typically presented 
in terms of the asymptotic behaviour of the tail: 

Mr,t)~t- 2e - 3 , cj> i (u,v)~v- u - 3 , rMu,v)~u- 2 '- 2 (11) 
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for Schwarzschild, where 2M < r < oo is fixed in the first formula, u > v in the 
second, and v > 2u in the third, and 

<Mr,;)~e- rft , ^(u,«)~e-^, ^(t!,!;)^-* (12) 

for Schwarzschild-de Sitter and ^ > 1, where r& < r < r c is fixed in the first 
formula, and it > w in the second, and v > u in the third. 

At first glance, statements (|12p may appear stronger than what is actually 
proven in Theorem 11.21 As quantitative statements of decay, however, state- 
ments (HU) and US} are in fact much weaker than what has now been mathe- 
matically proven, here and in [12] . For, rewriting, in particular, the first formula 
of (JT2J) as 

\Mr,t)\ <C,(r)e~ cl \ (13) 

then there is no indication as to what Ci(r) depends on, indeed, if there is 
any bound on Ci provided by some norm of initial data, and if so, what is the 
behaviour as i — > oo. This does not concern a mathematical pathology, but 
is intimately connected with the physical effect caused by the photon sphere. 
Indeed, a geometric optics approximation shows easily that if ()13|) is to hold 
and if Ci is to depend, say, on the initial energy of the spherical harmonic, then 
Ce — > oo as I — > oo. It is the rate of this divergence that would then determine 
the decay rate (if any) for <f>. 

If one is interested in quantitative statements of decay, a statement like (|12l) 
provides no more information than the statement 

lim (j>t(u, v) = 0. (14) 

(u : v) — >(oo,oo) 

It is worth noting that the above statement at the level of individual spherical 
harmonics, together with the (uniform) boundednes^l result 

H<CsupM+CEj(^9j), (15) 
can indeed be used to show, for fixed r, the statement 

lim U-<t> o )(u,v)=0, (16) 

(u,v) — >(oo,oo) 

for the total </>. This can be termed the statement of (uniform) decay without a 
rate. 

Thus it is truly only (|16p . and not the results of [12] or Theorem II. H that 
can be said to be suggested by heuristic and numerical studies. 

Results like (fTB"|) or even just (fT5|) are sometimes referred to as "linear sta- 
bility" in the physics literature^ One should keep in mind, however, that were 

6 The result (1 1 5 I t was shown for Schwarzschild in fundamental work of Kay and Wald | 16| . 
Our 1121 gives an alternative proof not relying on the discrete symmetries of the maximal de- 
velopment. For Schwarzschild-de Sitter, the statement (I15H of course follows from Theorcin ll. il 
for any s > 1. We have not found another statement of this in the literature. 

7 Sometimes, even the statement Vrj < r < r c , limt^oo <f>(r, t) = is termed "linear sta- 
bility". Such a result does not even imply 11151 . It is in fact entirely consistent with the 
statement sup rg(r()ir . c)i6[0iOo) \<p(r,t)\ = oo! 
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(fTo) the sharp decay result, it would in fact suggest instability for Schwarzschild 
or Schwarzschild-de Sitter once one passes to the next order in perturbation the- 
ory. At very least, it would exclude all known techniques for proving non-linear 
stability for supercritical non- linear wave equations like ([TJ. It is only quan- 
titative uniform decay bounds with decay rate sufficiently fast, such 
that moreover the bound depends only on a suitable norm of initial 
data, which indeed can the thought of as suggestive of non-linear sta- 
bility. One should thus be careful in associating the heuristic and numerical 
tradition exemplified in [19] with the conjecture that black holes are stable. 

1.5 Note added 

While the final version of this manuscript was being prepared, an interesting 
preprint [5] appeared addressing a special case of the problem under considera- 
tion here with the methods of time-independent scattering theory. The special 
case where <f> is not supported at H + DH^ and H. + DO. is considered and quan- 
titative exponential decay bounds are proven for JL, T (1 JT il n v in the coordinate 
t, where one must restrict to E' = {t} x [r ,i?o]> f° r r b < Tq, R < r c . The 
bounds lose only an e of an angular derivative, but depend on ro, Ro, and the 
initial support of <f> in an unspecified way. The work [5] depends in an essential 
way on a previous detailed analysis of Sa Barreto and Zworski [3T] concerning 
resonances of an associated elliptic problem. 

For the special case of the data considered in [5], given that result, then the 
estimates of the present paper, in particular, those provided by the currents 
J? , JK , can be applied a posteriori to obtain uniform (i.e. holding up to the 
horizons) exponential decay bounds. 
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2 The Schwarzschild de- Sitter metric in coordi- 
nates 

We refer the reader to the references [5J [T31 [TS] for detailed discussions of the 
geometry of Schwarzschild-de Sitter. 

2.1 Schwarzschild coordinates (r, t) 

We recall that so-called Schwarzschild coordinates (r, t) map T>° onto (Vfc,r c ) x 
(—00,00), in which the metric takes the form ([3]). Let the choice of the t 
coordinate be fixed. Here < r\, < r c denote the two positive roots of the 
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equation 

2M A , 
1 -r 2 = 0. 

r 6 

The function r can be given a geometric interpretation 
r(p) = yj Area(^" _1 (7r(p)))/47r, 



(17) 



(18) 



where here 7r : M. — > Q is the natural projection. Thus r can be defined as 
a smooth function on all of M.. It is known as the area-radius function. This 
function also clearly descends to Q. 

The (r, t) coordinates degenerate along the horizons 7i + U7i~ and H + WH , 
on which r ~ rb, r — r c , respectively. 

It is immediate from the explicit form of the metric that the vector field Jj 
is Killing in T>° . This extends to a globally defined Killing field T on (Ai,g), 

which is null along TL + U Hr and 7i + U 7i , and vanishes along 7i + n 7Y~ and 



2.2 Regge- Wheeler coordinates (r*,i) 

We now proceed to define two related coordinate systems on 2?°. Let us denote 
the unique negative root of (| 1 T[) as r_ , and let us set 

_ d f 2M 
dr \ r 3 

and similarly k c , k_. Wc now set 



2k c 



r 

1 



+ 2^7 log 



Jl-1 



2« b 



- C* 



r 

1 

r6 



where C* is a constant wc may choose arbitrarily. For convenience, let us 
choose C* so that r* = when r = 3M, the so-called photon sphere. We call 
the coordinates (r* , t) so-defined Regge- Wheeler coordinates. 



2.3 Eddington-Finkelstein coordinates (u,v) 

From Regge- Wheeler coordinates (r*,i), we can define now retarded and ad- 
vanced Eddington-Finkelstein coordinates u and v, respectively, by 

t = v + u 

and 

r* = v — u. 

These coordinates turn out to be null: Setting \i = — + ^Ar 2 , the metric takes 
the form 

—4(1 — /i)dudv + r 2 da§2. 
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We shall move freely between the two coordinate systems (r*,t) and (u, v) 
in this paper. Note that in either, region T>° is covered by (—00, 00) x (—00, 00). 

By appropriately rescaling u and v to have finite range, one can construct 
coordinates which are in fact regular on 7i ± and 7Y ± . By a slight abuse of 
language, one can parametrize the future and past horizons in our present (u, v) 
coordinate systems as 

^ + = {(oo^Ok.el-ocoo), ( 19 ) 

W = {(w,-Oo)}«£(-ao,oo]) (2°) 

H + = {{u, oo)} ue[ _ 00 , 00 ), (21) 

T-L + = {(-oo,u)}„ e[ _ 00i0o) . (22) 

Under these conventions, the statements of Theorem 11.11 can be applied up to 
the boundary of T>. 

2.4 Useful formulae 

Finally, we collect various formulas for future reference: 

2M 1 , 2 

M = + o Ar 2 , 

r 3 

9u V = (g uv r 1 = -2(1 -m), 
9„r=(l-/i), 9 u r=-(l-/i) 
dt = dv + du, dr* = dv — du, 

dt 2 9m 

d _l f d d 
dr* 2 \<% 9it 

dVol_\4 = 2r 2 (1 — /i) c?u cM§2 , 



<iVoZ t=CO n S t = r 2 y/l - fi dr* cL4§2, 
Oyj = V Q V Q ^ = -(1 - v)- 1 (a t V - r- 2 d r ,\r 2 d r ^)) + f A f A ijj. 
Here J7 denotes the induced covariant derivative on the group orbit spheres. 
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3 The energy momentum tensor and compatible 
currents 

As discussed in the introduction, the results of this paper will rely on L 2 -bascd 
estimates. Such estimates arise naturally in view of the Lagrangian structure of 
the wave equation. We review briefly here. 

Let <f> be a solution of D s (/> = 0. In general coordinates, the energy-momentum 
tensor T a p for <fi is defined by the expression 

T a p{4>) = d a 4>df3(j)- ^g aj3 g' 1& d 1 (t)ds4>. 

The tensor T a p is symmetric and divergence-free, i.e. we have 

V Q T Q/3 = 0. (23) 

For the null coordinate system u,v,x A 7 x B we have defined, where x A , x B 
denote coordinates on S 2 , we compute the components 

T vv = {d v (j)) 2 , 

Tuv = -\guv\f4>\ 2 = (i-m)IWI 2 - 

Here the notation iy-01 2 = 9 AB dA^pdBi/j — r ~ 2 \d^\ 2 dcr - Note moreover that 

IWf = r- 2 Ef=il^Vf- 

Let V a denote an arbitrary vector field. Let Tiy 13 denote the deformation 
tensor of V, i.e., 

nf = ^(V a V p + V' a V a ). (24) 
In local coordinates we have the following expression: 

T a p(<j>)T? = J^ZTJa ((du4>) 2 d v {V v {l - ^r 1 ) + {dA) 2 d u {V u (l - m)- 1 ) 

+\f4>\ 2 {d u v v + d v v u )) - ^-(K - v v ){\f<t>\ 2 - 4> a (t> a )- 

' 2r 

Set 

Jl = T a0 V a . (25) 

The relations (23]) and ((24]) give 

and the divergence theorem applied to an arbitrary region 1Z gives (fTOl) . 

Identity (JTOJ) is particularly useful when the vector field V is Killing, for 
instance the vector field T defined previously. For then, K T = and one 
obtains a conservation law for the boundary integrals. Moreover, when dlZ 
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corresponds to two homologous timelike hypersurfaces, the integrands Jj(^>)n /i 
on the right hand side of (|26[) when properly oriented are positive semi-definite 
in the derivatives of <f>. 

Were the vector field T timelike in all of V, then by applying (fT0|) to </>, 
fli4>, etc., one could show the uniform boundedness of all derivatives of 4>. Since 
T becomes null on Ti + U Ti , the integrand does not control all quantities on 
the horizon. It is for this reason that even proving uniform boundedness for 
solutions of Ogtp = on V is non-trivial. (See [16].) 

For <f> a solution to O g cj) = 0, the 1-form (</>) defined above has the property 
that both it and its divergence V M depend only on the 1-jet of <fi. Following 
Christodoulou [5], we shall call one- forms J M and their divergences K = V^J^ 
with the aforementioned property (thought of as form-valued and scalar valued 
maps on the bundle of 1-jcts, respectively) compatible currents. 

4 Constants and cutoffs 

4.1 The special values r i5 Ri 

In the course of this proof we shall require special values r», R%, i = 0,1,2, 
satisfying 

-oo < ^r* < 2r\ < ^r{ < 2r* 2 < < 2R* 2 < -RI < 2R\ < ^R < oo. 

Eventually, specific choices of these constants will be made, and these choices 
will depend only on M, A. Constants ri, Ri are in fact only constrained by 
Lemma 17.3.11 Constants t\ , R\ are constrained by the necessity of satisfying 
Proposition 18.3. ll of Section HOI 

To choose tq, Rq, on the other hand, is more subtle, as we will have to keep 
track of a certain competition of constants as ro, Rq vary. We adopt, thus, the 
convention described in the next section. 

4.2 Dependence of constants C, E, and e on r i5 Ri 

In all formulas that follow in this paper, constants which can be chosen inde- 
pendently of 7g, Rq shall be denoted by C. Constants C will thus depend on 
M, A, rj, Ri, for i = 1,2, and, after r^, Ri have been chosen, will depend only 
on M, A. 

Constants which depend on M, A, 7q and Rq and tend to as 7q — > — oo, 
Rq — > oo will be denoted by e. Finally, all other constants depending on M, 
A, Tq and Rq, will be denoted by E. Constants denoted by E in principle 
diverge as Tq — > oo, Rq — > oo. 

We will also use the convention A w B whenever C~ 1 A < B < CA with a 
constant C understood as above. 

In view of our above conventions, note finally the obvious algebra of con- 
stants: C ± C = C, eC = e, CE = E, eE = E, etc. 



18 



4.3 Cutoffs 

Associated to these special values of r, we will define a number of cutoff func- 
tions. It is convenient to introduce also the notation 

4.3.1 The cutoffs rj, 

Let rj : [0, oo) — > M. be a nonncgative smooth cutoff function which is equal to 1 
in [0,1] and outside [0,2]. 
For i = — 1, 0, 1, 2, define 

r]i(r*) = v( r * / r i) f° r r* < 
= v(r*/Rt) for r* > 0. 

Clearly rji has the property that r\i = 1 in [r* , R*] and r\i = in (— oo, ^r*_ x ] U 
(^R*-!,oo). Moreover, sup??- — ► as r* — > -oo, R* — ► oo. 

4.3.2 The cutoffs Xi an d Xi 

Now let x '■ (—00,00) — > K, x : (00,00) — > R be nonnegative smooth cutoff 
functions such that x is 1 in (— 00,— 1] and in [— 5,00), and x is 1 in [l,oo) 
and in (—00, i]. 

Now for i = —1, 0, 1 define 

Xi = X(~r*/r*) 
X t = X(r*/R*). 

Clearly, Xi nas the property that Xi = 1 i n ( _ °°i r *]7 an(I i n [2^*+i,oo). 
Similarly Xi has the property that Xj = 1 in [R*,oo), and in (— 00, 2R* +1 }. 

5 The hypersurfaces £ f and the region H(ti,t2) 

Let n, i?i be the special values announced in Section I4TT1 For all t, define 
E t = {t} x [n, Rt) U{(t-R{)/2} x [(t+i^)/2, 00] U [(t - rJ)/2, 00] x {(f + r£)/2} 
and for any £2 > ii, define 

ftfa.ta) = J+(s tl )n J-(E t2 ). 
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The diagram below may be helpful: 




We shall repeatedly apply the identity (fTU)) in the region TZ(t\ , t 2 ). We record 
below its explicit form: 

= / K{4>) + f J^K - / -WK (26) 

+ I JM n>1 + [_ JM n,x - 

JH+n{vi<v<v2} JH n{«i<«<« 2 } 

Here, n = (1 — n)~iT whenever (u,v) belongs to the space-like portion of St, 
and the measure of integration, call it dm, is there understood to be the induced 
volume form. Let us moreover define n = ^ and n = -J^ whenever (u, v) belongs 
to the v = const and u = const portions of S t respectively. With this choice, 
the measure of integration dm in the respective null segments is understood to 
be given by 

dm v=const = r 2 dA S 2du, dm u=const = r 2 dA&dv. 

All integrals over E 4 that appear in the sections that follow are understood 
to be with respect to the measure dm defined above. The integral over TZ(ti,t 2 ) 
is to be understood to be with respect to the spacetime volume form. Sec 
Section El 

6 The main estimates 

In this section we will give a geometric statement of the main estimates. 
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6.1 The vectorfields N, N and P 

Recall the cutoffs functions rji, XXi Xi; from Section [4731 Define the vector fields 



N = t^^- + t^^+T, (27) 
1 — fi ou 1 — \i av 

Pi = ^(l-^-Va^., ( 28 ) 
N=(r- 3M) 2 N. (29) 

For convenience, let us denote P 2 by P. 

The above vector fields will provide the fundamental directions in which the 
energy momentum tensor T M „ is to be contracted and/or </> is to be differentiated 
in the definition of the fundamental quantities appearing in the main estimates. 
See Section [6721 below. 

The coordinate-dependent definitions given above notwithstanding, the im- 
portant features of these vector fields can be understood geometrically. All three 
are invariant with respect to the action of ^>t, the one-parameter group of dif- 
fercntiable maps V — > V generated by the Killing field T, N is future-directed 
timclike on S t , N is future-directed timelikc everywhere except r — 3M, where 
it vanishes quadratically, and P is supported away from the horizon and orthog- 
onal to T. 

6.2 The quantities Z^' p , Z N and Q 

Let T^y be the energy- momentum tensor defined in Section [3l Define the quan- 
tities 

< ,p ® = I + (m 2 ) , (30) 

J St 

z?(<)= / r„„(0iw, (31) 

Q4>(ti,t 2 )= [ 2 Z^ P (t)dt. (32) 



Jti 

The quantity is equivalent to the spacetime integral of the density q(4>) 
defined by 

q{4>) = (t 11u {4>)^^— + (P4>? 



in the sense of the formula 



Q*(fi,t 3 )« / 9(0), (33) 



understood with the conventions of Sect ion |4~2I We will make use of this equiva- 
lence often in what follows. Recall also that the spacetime integral on the right 
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hand side of ([55]) is to be understood with respect to the volume form. Sec 
Section [2^41 

Note that the quantity Z^ has integrand positive definite in d<f>. It is in 
fact precisely the flux through S t of the current (</>). The quantity Z^' P 
differs from Z^ in that control of the angular and ^-derivatives degenerates 
quadratically at r = 3M . Similarly, the integrand of Q</,(ii, t 2 ) also degenerates 
at r = 3AI. This hypersurface r = 3M is the so-called photon sphere discussed 
already in the introduction. 

6.3 Statement of the estimates 

The main estimates of this paper are contained in the following 

Theorem 6.1. There exists a constant C depending only on M, A such that 
for all t2 > t\ and all sufficiently regular solutions <fi of ^ g 4> — in TZ{t\, t 2 ) we 
have 

Q*(ti,ta) <CZ$(h), (34) 

3 

z£(t a ) < C Zf p {t 2 ) + C(t 2 - h)- 1 (z$(ti) + Q^ii, t 2 ) + Y, Qn**(fi, t a )) , 

i=l 

(35) 

Z$ r (t 2 )<C(Z^(t 1 ) + Q (t 1 ,t 2 )). (36) 
More generally than §36}, ifT,' C 7?-(ti , ^2) * s achronal then 

f T^Wn" <C(Z^(t 1 ) + Q <p (t 1 ,t 2 )). (37) 

These estimates will be used in Section [TT] to prove Theorems 11.11 and 11.21 
As described in the introduction, the proof of Theorem 16.11 will be accom- 
plished in Section [10] with the help of so called energy currents J M associated to 
the vector fields Xz, Y, Y and 9. We turn in the next sections to the definition 
of these currents. 

6.4 Discussion 

Were it Z N,P on the right hand side of (|34[) . or alternatively, were Q defined 
as the time- integral of Z N , then inequality (|34| would immediately lead to 
exponential decay in t for Q(i, t*) (cf. Lemma [TTTTTTJ . 

The appearance of Qn^ on the right hand side of (|35|) signifies that the esti- 
mates "lose" an angular derivative. At the level of any fixed spherical harmonic 
(j>l, estimates (j34|) and ([35]) lead immediately to exponential decay for Q<^(i,i*) 
and Z^\ The nature of the loss of angular derivative in (|35|) means that for the 
total Q0(t, t*) and Z^ , one can only obtain polynomial decay in t, where the 
bound on the decay rate exponent is linear in the angular derivatives lossed. 
Exponential decay for <j> would be retrieved if the "loss in angular derivatives" 
in estimate (|35p were logarithmic. See the dependence in I in Lemma 111. 1.11 
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7 The J x family of currents 

Wc define in this section a family of currents, all loosely based on vector fields 
parallel to The role of these currents in capturing the role of the "photon 
sphere" has already been discussed in the introduction. 

7.1 Template currents jj* for an arbitrary V = f-Jjp; 
Let / be a function of r* and consider a vector field 

_d_ 

dr 

Define the currents 

JVW = V(^ + J(/ , + 2 1 7 ^/)a A1 (0) 2 -ia Al (/' + 2i-^/) 



y=m-!L- (38) 



l d lf + 2 1 -Jlf 



1 , 1 



2^-2— r^^V. 



^• 4 (0) = T^WV + ^'d^tf-^f'tf, 



and the divergences 
Note the identities 



/x' 1 - ^ r-3¥ 



2(1 - M ) r 
1 / //' aA _ M ( 3 8M\ . MA 2A 2 r 



(39) 



2r \ 1 — A* r ) r 4 \ r J 3r 2 
We compute 



^2/' + 4——^-/^ (f) a (f) a , (40) 



4 V r 
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r 4 

-i(r^ / " + ' / ")'- 

The expression J^' 3 is not a compatible current in the sense of Section [3l 
since K v,s depends on the 2-jet of </>, but it can be treated as such when restricted 
to eigenfunctions of /f\. 

7.2 Discussion 

The relation of the photon sphere to currents based on vector fields V of the form 
f(r*)d r * is most clear upon examining the modified current J^' 1 and noting the 
coefficient of \y<j)\ 2 in K Vtl vanishes precisely at r = 3M in view of ([39]) . This 
indicates that if one is to have say K V ' 1 > 0, the function / must change sign 
at r = 3M, and the control of the angular derivatives must degenerate at least 
quadratically. 

The task of choosing a suitable / is simplified by passing to the further 
modified current J V ' 2 which effectively "borrows" positivity from the d r -4> term. 



i\l-n J r r(l-n) J (l-//)r 



fj,'(l - pi) 



4 1 - /t" 



+ / 



M(. _8M\ MA_2AV\ . 2 



r J 3r 2 



(l-M) 



41-/x- 
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Finally, one can take advantage of the further flexibility provided by chosing 
separately V for each spherical harmonic, where now, after passing to the current 
J^' 3 , the O'th order terms are united with the angular derivative terms. In 
Section [7.31 we shall construct a current 

We do not in fact ensure that K X ' 3 > everywhere, but rather, only in the 
region ro < r < Rq. The region near the horizon will be handled with the 
help of the spacetime integral terms controlled by the currents J Y , J Y to be 
discussed in the next section. 

Once an initial positive definite spacetime integral (albeit modulo an error) 
is constructed, other quantities can be controlled with the help of auxiliary 
currents. These are defined in Section 17.41 It is there where we also use the 
current template J V ' 4 . 



7.3 The vector fields X e 
7.3.1 The case £ = 

Define 

fa = -r- 2 , 

and set X = fod r * . Let <f>o denote the O'th spherical harmonic of a solution 
U g (f) = of the wave equation. Consider J*°' (4>o), K x °'°((j) ) as defined above. 
We compute 

K*°.°(0 O ) = 2r" 3 {9 ^ 2 . (41) 

Note that 

K Xafl (fo) > 0. 



7.3.2 The case I > 1 

Consider for now the expression K V ' 3 for a general vector field V of the form 
(|38[) , applied to a spherical harmonic 4>£ with spherical harmonic number i > 1. 
We have 

K V ' 3 (M - (r^( 9 -(^)) 2 -^T^« (42) 



+ 

Define 



/„, l-3M/r M( 8M\ MA 2A : 



Mr)=^+D 1 -y r + ^f3- 8M ^ ■ MA 



The following lemma is proven in Appendix |A] 
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Lemma 7.3.1. For all I > 1, i/iere exists a unique zero rh e of the function 
hg{r) in [rb,r c ], and there exist constants r\ < 0, R% > 0, depending only on M 
and A, such that 

r2 < r ht < R 2 . 

Moreover, lim^oo r/ l£ — ► 3M. 

Let r% , i?2 «ow 6e fixed, chosen according to the above lemma. 

Were the middle term on the right hand side of (|4"2")l absent, then we would 
have K v > 3 (<f>i) > for any / such that /' > 0, f(r h( ) = 0. 

The middle term of (|42|) vanishes if /"' = 0, but, with the requirement that 
f( r hi) = 0, in this case the function / cannot be bounded. By suitably cutting 
off a function linear in r* , we can ensure that K V ' z {cj){) > in r$ < r < Rq. 

Let rjo be the cutoff of Section 14.3.21 Define 




and the vector field Xg by 




(43) 



We have 



K Xe ' 3 (M > o 

in the region Tq < r < Rq. Moreover, in this region we have in fact, 



(44) 



^^ + (r-r h ^j<CK x ^ e ). 




(45) 



(Recall our conventions for constants C from Section [472] ) 



7.3.3 The currents J X ' 1 



Finally, for i — 1,2,3 define the "total" currents 



and their divergences 



K x >\4>) 



K x °>°(<t> ) +Y,K x ^{4>i) = VJ^*(0). 
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7.3.4 Controlling the error 

Besides obtaining nonnegativity for K X ' 3 (<fi) in the region ro < r < Ro, we need 
to understand the error in the region r < ?'o and r > Rq. It turns out that 
this error can be controlled bj@ tq{4>), where q{4>) is a slightly stronger quantity 
than the energy density q(<p). 
Define the quantity 

m = (^W^y^ + T] - 1 ^ 1 _ + ^ l) \r*\- 5 - X \n\ 2 +r)i{P<f>)*) (46) 



Here 7/_i, xi, Xi are the cut-off functions defined in Section 14731 Note that 

+ (r - Ulf ((d«fl 2 + IWI 2 ) + (3,.-0) 2 . 

Compare with 



(1-M) a Xl (l-A*) 2 

(r-3M) 2 ((a t 0) 2 + |^| 2 ), 



^)-X 1 7^+X l7 ^ + (9^) 2 (48) 



where q{4>) is the density of the main quantity defined in 
The inequality replacing (|44|) which holds globally is given by 

Lemma 7.3.2. The inequality 

K x >*{ct>) >-e [ q(<f>) (49) 
holds on all spheres of symmetry. 

Proof. It suffices to consider the regions r < ro and r > Ro- Relation (|47|) 
implies 

V-i(Xl +Xi)\ , 1 2 



Cg(^) > ( I + '-^^ ) Ir*!"*- 1 |^| 2 , 



in these two regions, while 

J§ 2 7^1 Js 2 



> 



The function is supported in the region [2rg,ro] U [Rq,2Rq] and obeys the 
pointwisc bound 

|/f(r*)|<qr*|- 2 , (50) 



8 Recall the conventions regarding constants e in Section 14.21 
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under our conventions for the constants denoted C. Note that rj-\{xx +Xi) = 1 
in the support of /"'. Therefore, 

/ K x *{<j>)>-C [ J—lr*^^-^ 
Js, 2 Js 2 1 — M 

holds on all spheres of symmetry. We obtain (|4"5)) with 

e = C(\r*\- 1+5 + \BZ\- 1+s ). 

□ 

7.4 Auxiliary currents 

We will also need several "auxiliary" currents. 

7.4.1 Auxiliary positive definite pointwise quantities 

Let us first define, however, certain auxiliary positive definite quantities. The 
auxiliary currents K will be seen to bound these quantities when integrated on 
spheres of symmetry in the region r < r < R . 
Define 



qi{<t>) 


= m(<t>) = miT^MNWii - ^y 1 + {p 2 4>f) 


qm 




q?'(<t>) 


= Vi(<P - M 2 \ 


q\{4>) 


= m(r-3M) 2 \f(f>\ 2 , 




= rn(r-3M) 2 (T<j>) 2 . 



Here rji are the cut-off functions defined in Section [431 Note that 

?i(0)«C(^(0)+^(0)+gf(0)), (51) 

qf+i < qf 

for x = 0, a, a', 6, d, and that 

r-H(£+l) [ (r-3M) 2 qf(fa)r 2 dA s i= [ q^r 2 dA 2 
Js 2 Js 2 

for all spheres of symmetry 

The currents to be described in what follows are motivated by the problem 
of bounding the positive definite quanitites whose latin superscripts they share. 
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7.4.2 The current J x ° 3 

Define 

fi{r*) = -\m{r*W-r hl )\ 

where 772 is as defined in Section 14.3.21 

Note that (//)"' = -1 on [r 2 ,R 2 ], ft(r ht ) = 0, and ff = for r* < 2r* 2 , 
r* > 2R*. 

Set X° = ffd r ' and define 

e>i 

The current J x ' 3 (0) will allow us to bound the spacctime integrals of the 
quantities an d 1o ( ( t ) )- See Lemma riO.4.11 At this point, we can see from 

(|42p and (|45[) that for each ^, the pointwise bound 



tfifa) + gS {4>i) + t{t + l)(r - r h< ) 2 gg (&) < CK^fa) + CK x ?< 3 (<f> e ) (52) 
holds in ro < r < Rq. 

7.4.3 The currents J x " < 2 , J x "'° 
Define 

f b (r*)=V2(r*)(r-3M) 
Set A 6 = and define as before the currents J x "' 2 ((f> - <f> ), J X "'°((f>) and 

The current ,2 (4>— 4>o) will allow us to bound-in addition to the previous- 
the spacetime integral of the quantity q\{4>)- See Lemma TlPA 21 At this point, 
we can deduce from (1411). (l42l and (15211 the bound 



(«8(0 + 9o'(0) + 9o(0)) r 2 dA S 2 (53) 
(CA X < 3 (0) + Ci( r ^) + CK xb ' 2 (q> - O )) r 2 dA s2 



< 

Js 2 

on each sphere of symmetry with r* £ [r$, Rq]. 

The current J x ' o (0), in conjunction also with J x ,o (0o) an d J x ' 4 (0 — 4>o) 
to be defined below, will be used to help bound the quantity Z N in Proposi- 
tion [TUXTJ For now, note that from (gDJ), |[4T|). |g2J) and (JS3J), the estimate 

(go a (0) + 9o'(0) + 9o(0) (54) 

- C( m {l - /z)(2(r - 3M)/r + 1) + r? 2 (r - 3M)) c^^Vjr 2 ^, 

< / fc^ 3 ( ( /)) + C'^ xa ' 3 ( ( /)) + C^ b ' 2 (0- ( / ) o) + C^' O (0))r 2 dA s2 
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holds on each sphere of symmetry with r* £ [tq , Rq] . Note that the left hand side 
of (I54p a priori does not have a sign, in view of the term containing d^(f>d^4>. Af- 
ter defining J' YC '°(0 O ) an d J xd ' 4 ((j) — </>o) below, we shall be able to improve |54|) 
with ([57]), 

7.4.4 The current J x<= 

Define 

/ c (r*)-r 2 . 

Set A c = / c g^, and define as before J* c <°, i^ c <°. 
We have that 

K xc >°(0 o ) = 2r(d t M 2 - (55) 

7.4.5 The current J xd ' A 

Let us then finally define 

fV)=m(r*)(r-3M) 3 . 

Set X d = f d TrT and define as before the currents J x ,4 and K x ' 4 . 

The currents J x '°(<fio) and J x ,4 (0 — <f>o) will allow us to bound the space- 
time integral of the quantity qf (<fi). See Lemma riO.4.31 For now, it follows from 
(|4"0|) , (|5"3")) and (|55|1 that on each sphere of symmetry with r* E [tq, Rq], we have 

J a + tf'fa) + 8§(0) + g?($) r 2 dA s , (56) 

< jf (CX X ' 3 (0) + CK x °- 3 (^) + CK xb ' 2 (q> - O ) + Cif xc '°(0 o ) 

+C/v xd ' 4 (0-0o))r 2 dA S 2. 

The currents J xb ' (4>), J x "'°((j)o) and J xd,Ac {4> ~ 4>o) together allow us to 
estimate, on each sphere of symmetry for r* G [rQ,i?g], the quantity 

/ C V2 (2(r - 3M)/r + l)((d t cb) 2 - \fcb\ 2 ) r 2 dA s . (57) 

< { (cK x ^{dp) + CK xa - 3 ((j)) + CK xb - 2 (<j)- <f> ) + CK xbfl {<j)) 
Js 2 v 

+CK XC '°( ( f )Q ) + CK xdfi {<j> - O )) r 2 dA S 2. 
Again, the left hand side of (|57| does not have a sign. 



30 



8 The currents J and J 

In this section, we define the currents and Jf[ , associated to two vector fields 
Y and Y , supported near the black hole and cosmological horizons, respectively. 
The role of these currents in capturing the "red-shift" effect has been discussed 
in the introduction. 



8.1 The vector fields Y and Y 

Let < 5 < 1 be a small number, and define functions a(r*), (3(r*) as follows. 
Recall the cutoff functions \i, V-i from Section |4. 3. 2i and define 

«(r*) = Xi(r*)(2 - /i + ^(r*)^*!-'), (58) 

(3(r*) = 2r- 2 (2M/r b 2 - 2Ar b /3)- 1 X i(r*)(l - fi + r,-i(r*)\rT S )- (59) 
We have a > 0, (3 > 0. To see the latter, note that 



dr 



\ r 2 3 



> 

r=r b 



since r_ < r b < r c are three non-degenerate roots of (1 — fi), which is non- 
negative on [rt, r c ]. The above expression approaches as M and A tend to the 
extremal values. 

Similarly, define functions a, (3 as follows. Recall the cutoff Xi from Sec- 
tion [43731 and define 

a(r*)=xi(r*)(2- f i + r,- 1 \r*\- 5 ), (60) 

P(r*) = 2r- 2 (2Ar c /3-2M/r c 2 )- 1 x 1 (r*)(l- V + V-i\r*\- S )- (61) 

Again, note that a > 0, (3 > 0. 
Set K to be the vector field 

1 — \i du dv 

and Y to be _ 

a(r*) d — d 

l — (i ov ou 

The application of cutoff ?y_i is in fact not essential for the arguments of the 
paper. The cutoff simply ensures that Y and Y have smooth extensions beyond 
H+ and U + . 
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8.2 Definition of the currents 

Define the currents 



and set 



We have 



jj(0) = V(0), 

jJ{<t>)=T^ v Y v {4>), 
K V (4>)=V»jJ(cf>). 



K Y (0) = T F W< 



(dut) 2 ( (2M 2Ar\ A 

m H~ " — J a ) + W^jf (62) 



2 \ 1 — /i 1 — /i y r V 1 — // 



and similarly 



2(1 - M ) 2 V W 2 3 J J 2(1 -nf K ' 

- ±\f<f>\ 2 (a' - (0(1 - n))') + 1 (a - /3(1 - //)) 3„<W. 

8.3 Discussion and the choice of r l5 i?! 

In this section we exhibit the "red-shift" property of the currents J Y , J Y . This 
is essentially contained in Proposition 18.3.11 below. In the context of proving 
this proposition, we will choose the constants r\, R\. 

Note that the polynomial powers in the definitions (|58|) (|59| would be un- 
necessary had they not been present in the definition ([46]) . (Their presence in 
(|46|) is in turn necessitated by our application of (|50| .) 

A slightly unpleasant feature of these polynomial decaying expressions is 
that, if left bare, they would lead to vector fields Y, Y which fail to be C 1 at 
TL + ', TL , respectively. This would not in fact pose a problem for the analysis 
here. We prefer, however, to introduce a cutoff in (|58|) - (j6Tj) to emphasize 
the geometric nature of all objects involved in the proof. 

Proposition 8.3.1. For r* sufficiently small and R\ sufficiently large, depend- 
ing only on M , A we have 

CK Y (4>) > q(cj>), CK Y (4>) > q(4>) 
in r < n and r > R±, respectively. 
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Proof. We give the proof only for K Y . The proof for K Y is similar. 
Recall the functions a, /3 from Section [51 Let us denote by 

76 = 2r b - 2 (2A//r b 2 - 2Ar 6 /3)" 1 . 

In the region r <r\, we have 

a = 2 - » + V ^(r*)\r*\- S , f3 = lb {\ -fi + 77-1 (r*) |r* | " 5 ) . 

As a consequence, for all r\ < 0, 

a ' = (7? - ^ Ar ) (! - + V-i\rT S + W-i \rT 5 , 
/OAF 9 \ 

P' =T>(yr- 3 ArJ (1 - M) + W_ikT 4 + 7&<fy-ikT 
Recall that 

/2M 2 , \ 

and note in addition that rj'_ 1 > in the region r < r\. 

Recall now the expression (|62p for K Y . In the region r < r\ , we compute 

. a + ,_ 1 |,r i ) f^-?^ -cm- 



l-<5 



\ r 2 3 / \ r 2 3 

a' (/3(l-/i))' . M , . ... „ /2M 2 



(/i + 7b(l - A*) - 76^?-i(r*)|r*|- d ) — - -Ar 
1 — yul — /i \r^3 

+ ((1 - /i)- 1 - r fb)(v'-i\ r *\~ 5 + Sv^lrT 1 - 5 ), 
1(t^-P) = r- 1 (l- M )- 1 (l + ?? _ 1 (r*)|r*|- 5 - 1 ) + l 

- lb (l-f,)- lbV ^(r*)\r*\- S ). 



Note that 

r _1 (l - /i)" 1 < -(2M/r 2 - 2Ar/3)(l - ^)~ 2 + r" 2 (2A//r 2 - 2Ar/3)~ 1 . (64) 

It follows that n can be chosen sufficiently close to r b , where the choice depends 
only on M, A, such that, in the region r < n, we have 

ck*w > (J^s + (1 + T^i^r 5 - 1 ) ((^) 2 + \n\ 2 )) ■ (65) 

In deriving (|65p . we have used the Cauchy-Schwarz inequality to bound the 
d u (f>d v (f> term. (The fact that the constants work out is ensured by the limiting 
inequality ([64)) . It is here that the presence of the constant factor 75 in the 
definition of (3 is paramount.) 
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It now follows from (j4"T)) that, in the region r < ri, 

CX y (0) >g(^), 

as desired. 

Note finally, that, in conformance with our conventions of Section 14.21 the 
constant C above does not depend on ro, Rq, despite the appearance of the 
cutoff rj-i. 

□ 

Henceforth, let n , R\ be chosen so that the conclusion of the above propo- 
sition holds. 

The special values r±, R±, r%, R2 now being fixed, constants denoted C now 
depend only on M, A. 

9 The current J e 

Let £ : [0, 1] — > [0, 1] be a cutoff function such that 

C(x) = 1, x > 3/4, 

CO) = 0, x < 1/2. 

Define 

C(t 1 ,t 2 )(r)=C((T-ti)/(ti-i2)). 
Let 8 be the Heaviside step-function and let to be the vector field 

6 = (fl(rI-r*)C( tll t a )(2«-rI)+fl(r*-rI)fl(i^-r*)C( tllta )(*) 
+ 0(r*-i?nC (tl , t2) (2u + i?D) T, 

and define as before the currents J e and if . Despite the appearance of the 
Heaviside function, is a C 0,1 vector field. It is of the form = £T, where £ is 
a spacctime cut-off function adapted to the C 0,1 foliation S t , constant on each 
leaf Si . We have 

= -^T3^) (pm - ^)C(t 1)t2) (2« - rn (W) 2 + i(i - n)\n\ 2 ^) 

+26(r* - r^W - r*)C (tl , ta) (t) {{d t <pf + (8 r ^) 2 + (1 - »)\fcj>\ 2 ) 
+%* - iZf)C| tl|ta) (2« + i^) ("(a^) 2 + 1(1 - f,)\n\ 2 
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10 Proof of the main estimates 

10.1 Auxiliary integral quantities 

Let us define the auxiliary quantities 



Q£*(*i,*a) = 


Jn(t ly t 2 ) 


Q<p(ti,t 2 ) = 


! m 


F*(*i,*a) = 


[ j t (</>k + 1 jJ(0K 

Jn+n{v 1 <v<v 2 } JH + n{u 1 <u<u 2 } 







where qf{4>) is as defined in Section [7.4.11 with x = 0,a, b, . . ., and q(4>) is as 
defined in Section [7. 3. 41 

10.2 Preliminary inequalities 
10.2.1 Auxiliary quantity inequalities 

For the boundary terms, note first that on the support of r)i, 

1 =T. 



On the other hand, in the expression (|27[) . each term gives a nonnegative con- 
tribution to , and we have on E t 

J?n» « Xl 1M! +Il iM! + ( dt ^ + ( dr ^f + (1 - ( 66 ) 
^ 1 — 1 — M 

As a consequence, 

Z^&^CZ^), (67) 

and 

Z£(ti)<CZ£fe). (68) 
10.2.2 Boundary term inequalities for currents J 

In this section we address questions of size of the boundary terms generated by 
the currents J x , J A "° , J x \ J x " , J xd from Section[7]and J y , J F from SectionH 

Proposition 10.2.1. For J = J X ' l {4>), J x °' 3 (0) ; J x ^ 2 {4> - (f>o), J xb '°{4>), 
J XC ' a ((f>o), J x ^ 4 (0-0o) we have 



<EZl(t), 
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/ _ ■//<£Fj(ti I i J ). 

JH+n{v 1 <v<v 2 } J'H + n{u 1 <u<u2} 

Proof. We shall here only prove the proposition in a representative case of the 
current 

j x >*{4>) = J Xo '°(<M + £j^ 2 W>a 

£>1 



defined in Section [7.3.31 
Recall that 



Since on the space-like part of St, we have n = (1 — fx) ? T, it follows that there 
T^)T^n v = -^=T tlv {<j>)T' i T v = - i (T„„(0) + 2T W (0) + T vv (4>)) , 

since in addition T = l/2(d/du + d/dv). On the other hand, since n = <9/<9u 
and n = d/dv on the null segments u=const and it=const of St, respectively, 

T^TW = T^($T" (J^j = ^(T uu (0) + T uv {<f>)), 

T„ v {4>)T»n v = (J^j = \{T VV {4>) + T U .M)- 

Since the space-like portion of T, t corresponds to r-values r± < r < i?i we have 
that 

Zj(*)« / {TuM + ^-Xi)Tuu(^) + (l-Xi)T vv m. 

Therefore, 

Z?(<)»/ ((l-Xi)(M) 2 + (l-Xi)(M 2 + (l-M)l^| a )- (69) 



On the other hand, 



JH+n{ui<ii<B2} Jw n{ui<ti<« 2 } 



{d^f + l -f_ (d u ^) 2 . 

H+n{v!<v<v 2 } 1 JH + n{u 1 <u<u 2 } 



2 



We compare now ([69]) with the expression for J s J^ ,2 n M . Start with the 
current 



3G 

Since d/dr* = 1/2(8 /dv — d/du), we infer that on S t 

f |j, Xo '°(0oK| < c[ {(i-xi)(duM 2 + (i-xi)(d v M 2 + (i-^)\no\ 2 ) 

For £ > 1, consider now 

J?" 3 fa) = -fiT^&i) [dr*Y-\ (fl + Z^ft) d^l+\d, U + 2^//) 4>l 

Recall that the functions f'^r*) = 1 on the interval [rg,-Ro] an< ^ vanish for 
r* < 2rg and r* > 2Rq. On the space-like portion of £ t , i.e., for r G [n,i?i], 
we have 

|J^ 2 (^K| < C(T uv (<t> £ ) +T uu (fa) +T vv (fa) +$) , 
while on the null segments of S t , we have 

< E (T uv (<p e ) + T uu (fa) + [dutpi) 2 + (1 - n)<$) 
for v = const, and 

< E (T uv (<j, e ) + T vv (&) + {du&Y + (1 - n)<ft) 
for u = const. Therefore. 

f |J** 2 (&)n"| < E [ (T„„(^) + (l~Xi)Tuu(fa) + (1 - Xi)T vv (<t>i) + (1 - Mt) 

<E f ((l-x 1 )(du^) 2 + (l-xi)(d v 4>e) 2 + (l-^(\fM 2 + 4))- 
Summing over I and using the identity 

e(e + i) [ 4^4= / \fM 2 r 2 dAl 

we obtain 

/ I J^ 2 (0)n"| < E ( ((l-Xi)(9^) 2 + (l-Xi)(^) 2 + (l-M)|y0| 2 )<-EZj(t), 
as desired. On the other hand, 

/ jf 2 w=\ f^' (^o)+E^ i2 (^)] (Tf- 

J-H+n{ui<ti<D2} J«+n{» I <K» 2 } \ J \ ov J 
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On the horizon H + , we have (1 — /i) = 0, and thus T uv = 0, f' e = fl = and 

{iy=- ay = 

As a consequence, since < E, 

<E [ (O^) 2 <EF^(h,t 2 ), 

Jn+n{v 1 <v<v 2 } 

as desired. Similar arguments give the inequality on the horizon 7i + , as well as 
the inequalities for the other currents of the statement of the proposition. □ 

Proposition 10.2.2. For J = J Y ((f>), J Y '(<£), we have 

0< / J^<CZf p {t), 



JH+D{v 1 <v<v 2 } 



H+n{v 1 <v<v 2 } 



Moreoi 



H+n{v 1 <v<v 2 } 



J/fOK > o, 



H + n{vi<V<V2} 



jJ(0K = o, 



77 + n{«i<u<ti 2 } 



J Y {4>)n» > 0. 

'■H i "n{iti<u<u 2 } 

Proof. Once again, we shall here consider only the current J Y (<fi). The consid- 
erations for J Y are practically identical. 

By the construction in Section the support of J Y (4>) is contained in the 
region r* < \r\. This immediately implies that J y (0)|^+ = 0. Moreover Y 
is a bounded future-directed time- like vector field in the region r* < r* < |rf , 
which implies that there we have 

o< Jj(0K = r^(0)r^ 

< CT ia ,(<j>)T'*r 

= CJj(</>K < CJ*{<j>)n». 
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Restricted to the support of Y, the remaining part of St is contained in a null 
segment v = const. Thus, we have 

JJ(0K = (J^y = j^~T uu + I3T UV 

(a u 0) 2 + /3(i-/i)|^| 2 . 



f3 = 7b (l-» + r,- 1 (r*)\rT S ), 



The functions a, (3 are non- negative and in the region r < t\ are given by 

a = 2 - n + 77_i(r*)|r*|"' 5 , 
which implies that 

W) 2 



< 31 fa)n" < C 



1 — jj, 



(i-aOIW 



Comparing this to the expression for 3^ (<^)n^, given in (|66[) . we see that 

o< jJ(«^K<cj^(0K 

on the null portion v = const of S t . Since on this portion we have N = N , we 
finally obtain the desired inequality 

o< / jJ(0K<cz£' p (i). 

On the other hand, on 7i + wc have n — and 

' nil 



I -fx 



T UV +(3T VV = 2|^| 2 , 



since (3 = (1 - /x) = and T TO = (1 - yu)|W| 2 on H+. Thus 

/ JJ(0K=2/ |W| 2 . 

Jw+n{ui<t)<u2} Jw+n{Di<u<t>2} 



□ 



10.3 Applications of the integral identity for currents 

In this section we exploit the divergence theorem for compatible currents to 
relate various integral quantities. From Proposition 1 1 . 2 ."Tl Proposition 1 1 . 2 ."21 
and identity (|26[) . the following two propositions follow immediately: 



Proposition 10.3.1. For J = 3 X ' l {4>), 3 x ^ 3 ((j)), 3 X ^ 2 (4> - <j> Q ), 3 x ">°((j)), 
3 XC '°((t>o), 3 xd ' i {cp -4> ) andK = V^ 1 J M , we have 



K 



n(ti,t a ) 



<£(Zj(i 1 ) + Zj(t 2 )+Fj(t 1 ,t 2 )). 
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Proposition 10.3.2. We have 



I K Y {<j } )+ [ JJ(0K+ / JJ(0K < Czf p (tx), (70) 

Jn{t-L,t 2 ) JT, t2 JH+n{v 1 <v<v 2 } 

{ K Y (t)+[ JJ(0K+/ jf(ct>)n»<Czf p (h). (71) 

Jn(ti,t 2 ) J~s t2 JH + n{«i<u<ti 2 } 

Applying (|26p with the energy current gives us the following proposition 

Proposition 10.3.3. 

Z£(ia) + F£(f lf t a ) = Zj(ti) <CZ; (t x ), 
Proof. The statement follows immediately from A T = and (|67|) . ([68]) . □ 

Applying |26|) with the energy current Jj + J J + Jj , we obtain 
Proposition 10.3.4. Lei £' C 7?. 4ljt2 6e achronal. Then 



j; (0)n" < Z£ ( tl ) + C / —K (0) - A^ (0) 

£' J7i(4 1 ,t2)nJ-(S') 

In particular, 

Z$ (t 2 ) < Z$ (f x ) + C f -K Y (0) A F (0). 
Proof. The vector field T + Y + Y is timelike and one sees easily 



while certainly A T + A y + A Y = A y + A y . □ 

We have an alternative bound on Z^ (t 2 ) as follows 
Proposition 10.3.5. 

Z^(t 2 )<CZ^ p (t 2 )-C f K e (c/>). 

Proof. Recall that by the construction given in Section [9j the vector field 9 is 
future timelike, S|s t = T and 0|s tl = 0. The result now follows from the 
statement 

Z$(t 2 )<CZ^ p (t 2 )+C [ J^K 

Jz t2 

and the divergence theorem. □ 
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10.4 Bounding Q^, from 

In this section we establish the first key part of Theorem 16.1) that is to say, 
statement (|34|) . We begin with the following 

Proposition 10.4.1. 

Qi,*(ti,ta) < SZj r (t 1 ) + eQ^(ii,i 2 ) 
Proof. The proposition follows from the three Lemmas below: 
Lemma 10.4.1. 

QU + Q a i.4> < C I K x <\4>) + C [ K xa '\cp) + e Q . 
Lemma 10.4.2. 

QU<C f K X3 (q>)+C [ K xa < 3 (q>)+C [ K^^^-M+eQ^. 

•/*R(ti,t 2 ) JTZ(t 1 ,t 2 ) Jn(t 1 ,t 2 ) 

Lemma 10.4.3. 

QU < Of K x >\4>)+C( K x ^(J>)+c[ K x ">\<t>-ct>o) 
^*R(ti,t 2 ) JnittM) Jn(t u t 2 ) 

+ C f K xc >°(4>o) + C [ A- A " d - 4 (0-0 o ) + eQ0. 

Proof. The statements of the above three Lemmas follow directly from (|45|) . 
(j52")l . (|53"| and ([55]) , combined with the observation that, since the supports of the 
currents J x ,3 , J x ,2 and J x ,4 , as well as of the quantities qf (</>), are contained 
in the region {2r\ < r* < 2i?J}, since K x '°(4>o) is nonnegative, and since K X ' 3 
is nonnegative in {tq < r* < -Rq}, it suffices to apply Lemma T7. 3. 21 □ 

Proposition ll0.4.T1 now follows from Proposition ll0.3.T1 and Proposition ll0.3.3l 

□ 

Proposition 10.4.2. 

/ K Y {cp)+K Y {cp)<Czf P {t l ). 
Jn(t u t 2 ) 

Proof. This follows immediately by adding (f70|) . (fTTj) of Proposition 110.3.21 in 
view also of Proposition 110.2.21 

□ 

Proposition 10.4.3. 

Q*(ti,ta) <CQ h 4h,t 2 ) + C [ K Y (<j>)+K Y {<t>) 
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Proof. The statement follows from the inequality 

q(cj ) )<Cq 1 ^)+K Y ( ( f > )+K Y (cj ) ) 7 
which is a direct consequence of ([51]) , (|4"T|) and Proposition 18.3.11 □ 

It now immediately follows from Proposition 110.4.11 Proposition 110.4.21 and 
Proposition llO. 4/31 and our conventions regarding constants C, E and e that the 
following holds: 

Proposition 10.4.4. If is chosen sufficiently small, and R^ is chosen suffi- 
ciently large, depending only on M , A, then 

Qt(tx,t 2 )<EZ$(h). 

Henceforth, let ro , Rq be so chosen so that the conclusion of the above propo- 
sition holds. In particular, in what follows we shall need only make use of con- 
stants C depending only on M , A. 

Finally, since by flE}, (gHJ) 

Q*(*i,*2) < Q*(ti,ta), 
the statement (|34p of Theorem 16.11 follows immediately. 

10.5 Bounding from Z^' P , and Q n .^ 
Bound ([33)) of Theorem 16.11 is contained in the following 
Proposition 10.5.1. For all t\ < t 2 

3 

(t 2 ) < Czf p (t 2 ) + (ta - ti)- 1 c(zj r (h) + Q^h,ta) + ]T Qn,*(*i, *a)) . 

i=l 

Proof. By Proposition 110.3.51 

Z$(t 2 )<CZ*- p (t 2 )-C [ A' e (0). 

Jn{t lt t 2 ) 

The current K®((f>) was defined in Section [5J One easily sees 

< j^— fxi^ + + IW 2 + * (W) 2 + (3^) 2 )) . 

*2 - t\ \ 1 - /Lt 1 - /X / 

Comparing this with (|48|) , we infer that the statement of the proposition would 
follow from the estimate 

r 3 

/ m{{9t<i>) 2 + \?<i>\ 2 ) <c(z^(t 1 ) + Q^(t 1> t 2 )+5^Q n< ^(t 1> t a )). (72) 
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From (I4"51) we have that 



CQ^h,t 2 )> f (r - Ml) 2 \f4>\ 2 + (d r * 

JK(ti,t 2 ) 

3 



I(tl,t2) 



A one-dimensional Poincare inequality immediately implies that 
c(Q*(* 1) t 2 ) + y;Qn 4 *(*i,ta) J > / m\f4>?- 

\ i= l / M*l,t2) 

It thus follows from ([57)) that 

/ (cA x < 3 (0) + CK xa ' 3 (^) + CK xb > 2 (<j) - O ) + 



2 



> / miW) 



C f Q^(ti, t 2 ) + Qn«*(*i,*a)l • 



The desired (|72p now follows from Propositions 110.3.11 and 110.3.31 □ 

10.6 Bounding Z$(t 2 ) from Z^(ti) and Q 
The final statements of Theorem 16.11 follows from 
Proposition 10.6.1. Let £' C 7£(ii,t2) be achronal. Then 

J T^faJJW < C(Z% (fx) + Q (ti, t 3 )) . (73) 

In particular, 

Z$ r (t 2 )<c(z$ r (ti) + Q^(ti ) t 2 )). (74) 
Proof. By Proposition 110.3.^1 



j; < Z£ (t0 + C / - A y (0) - A (0) 

B' ^R(Si,ij)nJ-(E') 

Recall that the current K Y ((f>) (respectively, K Y (<f>)) is positive for r < ri 
(respectively, r > R\) and vanishes for r* > \r\ (respectively, r* < 
Moreover, comparing ([151) and ([6"2"[) easily implies the bound 



-K Y {<f)-K Y {4>) <Cq{4>) 

for rl < r* < ir*, and ^i?* < r* < \R\- The result now follows immediately. 

□ 
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11 Proof of Theorem 11.11 
11.1 Energy decay 

Proposition 11.1.1. There exist constants C , c, depending only on M , A, such 
that for all <f>£ solutions of D g (f>£ = in J + (Eo) n T> with spherical harmonic 
number I, then 

Z%(t) + Q^(tX)<CZl(0)e-^ 2 
for all t and all t* > t. 



Proof. This follows immediately from estimates (|34j) and (|35|) . in view of the 
following lemma, proved in Appendix 151 applied to the functions f(t) — Z^ t (t) 
a ndh(t)=Z$ e (t): 

Lemma 11.1.1. Let k, fco be positive constants and let f : [0, oo) — > R, 

g; [0, oo) — > M be nonnegative continuous functions satisfying 

h(t 2 ) + jf 3 f(r)dr < k (f(t 2 ) + (i a - ti)- 1 ^(<r) + ^ 2 jf 2 /(r)dr) 

/or aZ/ i 3 > £ 2 > ii > 0, and f < fco- Then there exists a constants c 
depending only on k, and a universal constant C such that 

hit) + / f(t)di < C(max{/i(0), ko})e~ ct/e2 



for all t and for all t* > t. 

□ 

Proposition 11.1.2. There exist constants C , c, depending only on M , A, such 
that for all <f>£ solutions of O g (f>£ = in J + (Eo) n T> with spherical harmonic 
number £ and for all achronal E' C T> PI J + (Eo), 

T^iWn" < CZ^0)(e- ta ^''/' ! +e- 2c »+( E ')/ f! ), 



in particular 



Z»(f)<CZj'(0)e- M / 



Proof. This follows immediately from Propositions 111. 1.11 and the inequality 
(EH) of Theorem O □ 



The energy decay statements of Theorem 11.21 now follow from the following 

Proposition 11.1.3. Let E be a Cauchy surface for M., and let <p, ip, ~Eo((p, ip) 
be as in the statement of Theorem [777] or Theorem Then, there exist con- 
stants C , to > depending only on M , A and the geometry o/E n J~(T>), such 
that 

E t0 c J+(s)nP, 
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holds. 



Proof. This is completely standard. We give a sketch to emphasize here too the 
role of compatible currents based on vector fields! Extend say Nq from Eo to 
an arbitrary future-timclikc vector field N in J~(£t ) n J + (£), and consider an 
arbitrary spacclikc foliation S T of this region by manifolds with boundary, such 
that 5_i = £, So — St, and dS T C J _ (£t) \ I~(St)- Consider the analogue 
of (26]) in J-(S T ) n J+(S_i). Set 



On the other hand, one easily sees that there exists a C depending on the 
geometry of our chosen foliation of the compact set J~(Et) n J + (E) such that 
for all r G [-1,0], 




In view of the fact that N is future timelike, we obtain 




K. 




It now follow that /(0) < e c /(-l). 

The result now follows by noting that 



f(-l)<CEo{<p,<p), 



Z?(fo) = /(0). 



□ 



11.2 Point wise decay 

The pointwise decay statements follow easily. 



A Proof of Lemma 17.3.11 



Consider the function 



A(r) 



= r 3 h i (r)=£(£+l)[l- — ) + 




3M 8M 2 MAr 2A 2 r 4 



3 9 



Recall that 
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We may thus rewrite the function as 



A(r) = i(l + 1) (l - ^Pj +^~+ 3AAfr - 2 M 2 . 



For I > 1, A is concave in 



r: 



The value of A at r c is given by 
A(r c ) = {e(l+l)-2)(l-—)- — +3AAIr c + 2-2^ 2 (r c ) 

= (t(£ +l)-2)(l-—)- — + 3AMr c , 

where we have used fJ-(r c ) = 1. The same formula evidently holds at r = with 
rf, replacing r c everywhere. 

Since r c > 3M, + 1) > 2 and Kr 2 c = 3 - we have 

, . . 3M„ o , 3M / 6M\ „ 

^ c > ( Ar c - 1) = 2 > 0. 

r c r c \ r c J 

From < 3, we obtain similarly 

Al . 3M / 6M\ 
A(r 6 )< 2 <0. 

By concavity the function A(r) then has exactly one zero on the interval (r^, r c ), 
and thus, so does hi(r) = r~ 3 A(r). 

The second statement of the lemma is clear from the final, which in turn 
follows immediately from the form of the function hi. 



B Proof of Lemma 111.1.1 



By replacing k with maxjfc, 1}, we may assume in what follows that k > 1. Let 
ti, . . . ,U be a sequence with 18k(£ 2 + 1) > tj — U-i > 9fc(^ 2 + 1). We can choose 

£ + 9fc(£ 2 + 1), i 4 + 18fc(^ 2 + 1)] 

such that by pigeonhole principle 

f(U+i) < k-\M 2 + 9)- 1 / f(r)dr. (75) 

Jti + (9£ 2 +9) 

Assumptions of the Lemma then imply that 

h{U) + £ +1 f(r)dr < kf(U) + ±h(ti-l) + I f f(r)d(r), 



4G 



/•t, + (18£ 2 + 18) i 

/ f{r)dr < kf(U) + -h(U-x) + £ / /(r)d(r). 

'ti + (9^ 2 +9) y 



Therefore, 



/(fi+i) < (9^ 2 + 9)- 1 f(U) + — M^-i) + 777: / f(r)d(T) 



9k v ' L ' 9k 



U-i 



Thus if 

/(*,•) < C2- 1 - 1 (76) 

fe(ti_i) + / ' /(r)dr < CA:2- 4+1 (77) 
•/ti-i 

then 

< C'(2~ t 3" 2 + 3~ 4 2- t+2 ) < C2- 1 - 1 

h(U)+J + /(r)dr < C{k2~ i - 1 + 3~ 2 fc2 _m ) < C , /c2" 4 . 

Now we have that (|76|) . ((77)) indeed hold for i = 1. where to = 0, and ti is 
defined so as to satisfy ([75]) , with C = max{/i(0), fco}- This proves that for the 
sequence ij, defined above, 

+ ^ /(r) dr < C k 2~ { < C k e " c V« 2 . 

Using the assumptions of the Lemma again we immediately obtain the desired 
statement. 
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